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ABSTRACT 



We investigate the molecular gas content and the excitation and fragmentation properties in the central region of the spiral galaxy 
Messier 81 in both the 12 CO(1-0) and 12 CO(2-l) transitions. We have recently observed the two transitions of CO in the M 81 center 
with A, B, and HERA receivers of the IRAM 30-m telescope. We find no CO emission in the inner ~ 300 pc and a weak molecular gas 
clump structure at a distance of around 460 pc from the nucleus. Observations of the first two CO transitions allowed us to compute the 
line ratio, and the average I21 //10 ratio is 0.68 for the M 8 1 center. This low value, atypical both of the galactic nuclei of spiral galaxies 
and of interacting systems, is probably associated to diffuse gas with molecular hydrogen density that is not high enough to excite 
the CO molecules. After analyzing the clumping properties of the molecular gas in detail, we identify very massive giant molecular 
associations (GMAs) in CO(2-l) emission with masses of ~ 10 s M Q and diameters of ~ 250 pc. The deduced N(H 2 )/Ico ratio for 
the individually resolved GMAs, assumed to be virialized, is a factor of ~ 15 higher than the standard Galactic value, showing - as 
suspected - that the X ratio departs significantly from the mean for galaxies with an unusual physics of the molecular gas. 

Key words. Galaxies: Interstellar Medium - Galaxies: Spiral - Interstellar Medium: Clouds - Interstellar Medium: Molecules 



1. Introduction 

The molecular gas content o f Messier 81, one of the near- 
est face-on spirals (SA(s)ab; Ide Vaucouletirs et"al ] 119911) . has 
always been a puzzle as seen in the first observ ations by 
ISolomon & de Zafral d 1975b and lCombes et al.1 (Il977l) . The CO 
emission appears very weak in this galaxy, especially in the 
central regions, and the molecular content seems to be con- 
fined to the HII reg ions in the spiral arms, between 4 and 7 
kpc from the center dSotomon et al1ll987HBrouillet et alj 1988, 



199ltlSage & W estpfahl 199 ll;ISakamoto et alj200ll) Two stud- 
ies dSolomon & de Zafral 1 1975c ICombes et al.l 1 19771) searched 
the molecular gas in the center of M 8 1 but only found upper 
limits or failed to detect any - probably due to the low CO emis- 
sion - so they did not draw any conclusions from this result. The 
fir st detection of the molec ular gas in the center of M 8 1 was 
bv lSage & WestpfahU (1199 ll) . who found t hat the peak intensities 
are a factor 4 lower than those observed by Brouillet et al. ( 1988) 
in the outer disk. After the weak d etection bv lSage & Weslpfahl 
dl99ll) . onlv lSakamotoetal.1 d200ll) have studied the CO emis- 
sion in the center of M 8 1 and also found a particular molecular 
gas physics. M 81 has also been studied in a ll components of 
its interstellar medium (ISM), from radio (e.g. lBeck et al.lll985t 

al. 



Bietenholz et all 1 1996b . th rough optical and UV (e.g. iHoet 
1996|), to X-ray bands (e.g. lPeUegrini et al.ll2000l) . 



In addition to the weak CO emission, the HI profile also re- 
flects a central deficiency in atomic gas. The HI gas is abun- 
dant only in the inner spiral arms of the galaxy (IVissedll980t 
lAllen et al1ll997l) . and the shapes and the motions observed in 
HI in the spiral arms quite clearly obey the predictions of the 
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density-wave theory. This strong density wave is attributed to 
the tidal interaction with the othe r galaxies of the M 8 1 group, 
especially M 82 and NGC 3077 dKaufman et al.lll989t) . 

In general, interacting galaxies like M 81 possess higher 
CO luminosity by alm ost an order of magni t ude than non- 
interacting galaxies dBraine & Combesl 119931 : ICombes et al.l 
U994tlCasasola et alJl2004l ). Then, M 81 seems an exception re- 
markably poor for its CO emission. For this reason, M 8 1 is a 
kind of prototype of CO-poor galaxies, as is Andromeda or the 
center of NGC 7331, a lesser extent. 

M 8 1 is also a good candidate for exploring the problem of 
the varying CO to H2 conversion ratio (X = N(H2)/Ico)- This ra- 
tio is well-known for varying subs tantially in dwarf galaxies, for 
which the metallicity is deficient (Rubi o et al.ll 1 993t iTaylor et al.l 
1998), while for galaxies like M 81 that are giants with normal 
metallicity, the problem of the low CO emission has not been yet 
explored much. The CO to H2 conversion ratio for M 8 1 is sus- 
pected - on the base of the observations that are already available 
- of departing signifi cantly from the me an (X = 2.3 x 10 20 mol 
ctrT 2 (K km s -1 ) -1 , Strong etaT 1988), because galaxies with 
low CO conten t, such as M 31, NGC 55, or LMC, show a much 
larger X factor (lNakai&Kunolll995l) . 

In this paper, we assume M 81 galaxy at a distance of 3.6 
Mpc and adopt H Q = 70 km s -1 Mpc" 1 . We present the CO(l-O) 
and CO(2-l) observations of the central ~ 1.5 kpc of M 81 car- 
ried out at the IRAM 30-m telescope. The primary goal of this 
work is to investigate the amount and distribution of the molec- 
ular gas component in the central regions. 

The structure of the paper is as follows. In Sect. 2 the obser- 
vations performed and the techniques of the data reduction used 
are described. In Sect. 3 we present the observational results, and 



2 



V. Casasola et al.: Molecular clouds in the center of M81 




Table 1. Fundamental parameters for M 81 (at a distance of 3.6 
Mpc adopted for M 81, 1" is -17.4 pc). 



R.A. 


Dec 


V 


V 


line 


(J2000.0) 


(J2000.0) 


(km/s) 


GHz 




09:55:33 


69:03:55 


-34 


115.271 
230.538 


12 CO(1-0) 
12 CO(2-l) 



in particular, the R21 = hi/hi li ne ratio is studied and analyzed 
for all the offsets observed that give information on the physi- 
cal properties of the gas in the M 8 1 center. The clumpy nature 
of the molecular gas, the virial equilibrium, and the X-ratio of 
the resolved individual clouds are discussed in Sect. 4. Our re- 
sults are compared with previous studies of the molecular gas in 
M 8 1 and discussed, giving their physical interpretation in Sect. 
5. Finally, in Sect. 6 we summarize the main results. 



2. Observations and data reduction 

The observations were made from 2 from 6 January 2006 and 
from 31 March to 1 April 2006 at the IRAM 30-m telescope at 
Pico Veleta, near Granada (Spain). We used two different config- 
urations for the two observing runs and we covered the central 
90" x 90", where at the distance of M 81, 1" = 17.4 pc. The 
receiver-cabin optic system of the IRAM 30-m telescope is op- 
timize d to observe with up to four different receivers simultane- 
ously dWildlll999l) . 

In the first observing run we used the A 100, A230, B100, 
and B230 heterodyne receivers of the telescope, which simulta- 
neously observed at both 1 15 GHz (A100 and B 100, at A 2.6 
mm, where A is the wavelength) and 230 GHz (A230 and B230, 
at A =! 1.3 mm). We used the 1 MHz back-ends with an effective 
total bandwidth of 512 MHz at 2.6 mm and an autocorrelator at 
1.3 mm. We also used the 4 MHz filterbanks with an effective 
total bandwidth of 1024 MHz at 1.3 mm. These arrangements 
provided a velocity coverage of 1300 km s _1 in both the CO(l- 
0) and CO(2-l) lines, and all the measurements were performed 
in "wobbler-switching" mode. This observing mode has a min- 
imum phase time for spectral line observations of 2 sec and a 
maximum beam throw of 240". The advantage of the "wobbler- 
switching" mode is to give very flat baselines without any ripple 
in most cases. The half power beam widths (HPBW) are 22" 
and 12" for the CO(1-0) and CO(2-l) lines, respectively, and the 
typical system temperatures were ~100 K at 115 GHz and ~400 
K at 230 GHz. The positions observed with A and B receivers 
were 81 for both lines. 

In the second obs erving run, we us ed the Heterodyne 
Receiver Array HERA dSchuster et alj|2004 . a focal array of 18 
SIS receivers, 9 for each polarization, tuned to the CO(2-l) line 
for M 8 1 . The 9 pixels are arranged in the form of a center- filled 
square and are separated by 24". The sampling was 6", and a 
homogeneous mapping procedure was used to regularly sweep a 
12x12 pixel map, filling the intrinsic square of 66" x 66". The 
typical system temperature was ~400 K. Also in this case we 
used the "wobbler-switching" observing mode, and the pointing 
accuracy was ~ 3" rms. The WILMA backend was used, provid- 
ing a 1 GHz bandwidth for each of the 18 receivers. The bands 
contain 2x465 MHz (=930 MHz) channels spaced by 2 MHz 
of resolution. The total bandwidth corresponds to 1300 km s _1 
with a velocity resolution of 2.6 km s _1 . The map realized with 
HERA included 144 positions. 



Fig. 1. The field of our observations (white box, 40 arcsec in 
size) superposed on a GALEX composite image of M 8 1 . In the 
figure, FUV is in the blue channel and NUV the yellow one. 



The line intensity scale throughout this paper is in units of 
T m b, the beam-average radiation temperature. The value of T m b 
is related to T* A , the equivalent antenna temperature -corrected 
for rear spillover and ohmic losses- reported above the atmo- 
sphere, by T] — T* A IT m b, where 77 is the telescope main-beam 
efficiency, 77 defined as the ratio between the beam efficiency 
B e ff and the forward efficiency F e ff, 77 = />',.,, //•',. . The beam 
efficiency is usually derived from continuum measurements of 
planets, while the forward efficiency is obtained by pointing the 
antenna at different elevations and measuring the T* A of the sky 
at each step. The sky temperatures are successively fitted by an 
exponential function of the air mass, and the forward efficiency 
is derived with the atmospheric opacity. At 115 GHz 77 = 0.79, 
while T} Beff = 0.57 at 230 GHz. 

The data were reduced with the GILDASQ software package. 
Some spectra with random, highly non-linear baselines were 
suppressed, and linear baselines were subtracted from all the re- 
maining spectra. In Table [TJ the fundamental parameters for the 
M 81 observations are summarized. 

3. Molecular gas emission 

The region observed with the A and B receivers of the IRAM 30- 
m telescope in the two CO lines covers the central ~ 1 .6 kpc of 
M 8 1 . The size and the location of the region covered by our 
CO observations are shown in Fig. Q] as a square superposed 
on a GALEX image of the galaxy. The CO emission, both in 
CO(1-0) and CO(2-l) lines, come s from different regions o f the 
galaxy nucleus. I n agre ement with Sage & Westpfahl (1 1 99 lb and 
ISakamoto et alj d2001l) . we find that the nucleus presents some 
regions devoid of molecular gas, but also that there are other re- 
gions in which the CO is clearly detected and the signal is strong 
(Fig.0. 

In particular, the central ~300 pc regio n is devoid of CO(l-O) 
emis sion in agreement, for instance, with iBrouillet et all (l988) 
and [Sakamoto et alj d2001l) . but there is a main molecular gas 
clump structure at a distance of around 460 pc from the nucleus 
in the northeast direction (Fig. [3]), which is very similar to the 

1 http://www.iram.fr/IRAMFR/GILDAS/ 
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Fig. 2. Map of the observed M 81 central region. Each spectrum has a velocity scale from -600 to 600 km s -1 , and a temperature 
scale from -20 to 70 mK. The positions are offsets relative to the M 81 nucleus assumed of coordinates RA/2000.0 = 09 /! 55 m 33 s , 
DEC/20000 = 69° 03' 55". Upper panel: Observations of the I2 CO(1-0) line. Bottom panel: Observations of the 12 CO(2-l) line 
convolved to the 12 CO(1-0) beam resolution (22"). 



maximum intensity detected bv lSakamoto et al.l d200ll) . In addi- 
tion to the CO arc in the northeast direction, we also detected 
some emission to the southwest one. 

3.1. HERA-receiver results 

The observations done with the HERA receiver cover the cen- 
tral ~ 1.3 kpc, similar at the central ~ 1.6 kpc reported above 



observed with the A and B receivers. These observations con- 
firm that the molecular gas emission in the M 8 1 center is at a 
low temperature. Fig. [4] shows the spectral map and Fig. [5] the 
integrated contours obtained with our observations. HERA ob- 
servations performed for the 12 CO(2-l) line also reveal "islands" 
of molecular gas in the nuclear region. There are regions com- 
pletely devoid of emission, but also central offsets with a faint 
but detectable emission. The pointings with signals > 2cr have 
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Table 2. CO(/=1-0) and CO(/=2-l) intensity lines where CO emission is detected (detection level > 2(f). 



Offset 
Aa (") 
(1) 


Offset 

A/3(") 
(2) 


7^(1-0) 
[mK] 

(3) 


Av(l-O) 
[km s~'] 
(4) 


7* 

J r m6 (i-o)dv 

[Kkms~'] 
(5) 


7^(2-1) 
[mK] 
(6) 


Av(2-1) 
[km s~'] 
(7) 


T 

J T mh (2 - l)dv 
[Kkms 1 ] 

(8) 


(9) 


-51.9 


-30.1 


29.04 


74.01 


2.15 


10.23 


83.21 


0.85 


0.40 


-47.2 


-19.1 


30.16 


59.01 


1.78 


12.73 


82.55 


1.05 


0.59 


-42.5 


-8.0 


28.34 


88.23 


2.50 


14.55 


115.60 


1.68 


0.67 


-37.8 


3.0 


16.97 


79.81 


1.35 










-36.2 


-23.8 


34.22 


59.89 


2.05 


16.67 


87.36 


1.46 


0.71 


-31.5 


-12.7 


38.18 


52.72 


2.01 


23.87 


73.12 


1.75 


0.87 


-26.8 


-1.7 


18.09 


150.30 


2.72 


20.39 


97.33 


1.98 


0.73 


-25.1 


-28.4 


31.67 


61.07 


1.93 


25.13 


52.02 


1.31 


0.68 


-20.4 


-17.4 


33.16 


47.82 


1.59 










-19.1 


47.2 


31.85 


56.16 


1.79 










-8.0 


42.5 


38.36 


96.15 


3.69 


16.16 


70.65 


1.14 


0.31 


-1.7 


26.8 


68.42 


102.50 


7.01 


40.16 


102.10 


4.10 


0.58 


3.0 


37.8 


27.99 


152.80 


4.28 


24.86 


116.50 


2.90 


0.68 


6.4 


-15.7 


12.52 


113.40 


1.42 










9.4 


22.1 


39.60 


78.91 


3.12 










14.1 


33.1 


34.47 


66.19 


2.28 


22.30 


80.29 


1.79 


0.78 


15.7 


6.4 


41.16 


44.16 


1.82 


18.33 


94.03 


1.72 


0.95 


22.1 


-9.4 


31.75 


96.58 


3.07 


20.61 


121.10 


2.50 


0.81 


26.8 


1.7 


21.40 


67.65 


1.45 


22.95 


58.55 


1.34 


0.93 


33.1 


-14.1 


21.76 


43.95 


0.96 










48.9 


-7.7 


58.60 


32.25 


1.89 


19.20 


52.02 


1.00 


0.53 


53.6 


3.3 


47.86 


55.29 


2.65 


35.44 


61.72 


2.19 


0.83 



The interpretation of the R21 = hi/ hi = J T,„b(2 - l)dv 
/ J T m h{\ - 0)dv line ratio is quite complex and usually simpli- 
fied with standard assumptions. The main assumption is to sep- 
arate the analysis of the R21 ratio into two limiting cases: the 
optically thin case and the optically thick case. The analysis of 
these two limiting cases for the CO molecule has shown two im- 
portant things. First, for high temperatures, the R21 ratio asymp- 
totically approaches a value of 4 in the optically thin case, while 
it only approaches 1 in the optically thick case. If 1 is the limit 
for an optically thick gas, R21 ratios > 1 are signatures of opti- 
cally thin CO emission. Second, in the optically thin limit, the 
R21 ratio does not drop below 1 until the excitation temperature 
drops below 8 K. This situation is improbable, but not impossi- 
ble, especially if the excitation of the gas is subthermal. 

The 12 CO lines are optically thick in most galaxies, includ- 
ing the Milky Way, but there are other galaxies where the prob- 
lem is still open, like M 82 -the inter acting companion o f M 81- 
(iLoiseauetalJI 19901) or NGC 3628 dReuteretal]ll99"ll) . When 
studying t he line ratio for a large s ample (118) of nearby spiral 
galaxies, Braine & Combes (1992) saw that subthermal excita- 
tion or different gas filling factors are only required for explain- 
ing line ratios less than 0.7. 

With line ratios of 0.7, 0.8, and 0.9, the excitation temper- 
atures are expected to be 7 K, 10.5 K, and 21 K, respectively. 
These low line ratios can be explained considering that at high 
gas densities in the heating mechanism the dust-gas collisions 
are more important than cosmic rays that ionize H2 or via for- 
mation of H2 on dust grains. When neglecting the heating from 
gravitational collapse or magnetic sources and assuming the 
gas in local thermal equilibrium (T ex - ThA iBraine & Co mbes 
(1992) derived a relationship between gas temperature and dust 
temperature at high densities. At gas densities of 10 4 cm -3 , 
which is required to approximately thermalize the J = 2 level 
of CO, and dust temperatures of 30 K, the equilibrium gas tem- 
perature is 10 - 13 K. Considering that in the clouds the dust is 
probably warmer in the vicinity of OB stars and lower elsewhere, 
the gas heated by the warm dust at 50 - 60 K will reach the ther- 




40 20 -20 -40 

RA offset (arcsec) 

Fig. 3. Integrated 12 CO( 1 -0) emission. The contour levels are be- 
tween 0.701 to 7.01 K km/s in steps of 0.701 K km/s. The beam 
of 22" is indicated at the bottom left. 



peak temperatures between 16 mK and 65 mK. There is good 
correspondence between the emission in CO(2-l) found with 
HERA receiver and with A and B receivers, as shown in Fig. 
[6] where the distributions found with the two observing configu- 
rations overlap. 

3.2. The CO(2-1) to CO(1-0) line ratio 

The combination of observations of the CO(1-0) and CO(2-l) 
lines allow us to compute the line ratio and to study the physical 
properties of the gas more in detail, such as excitation tempera- 
ture and optical depth. 
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Fig. 4. 12 CO(2- 1 ) map of the M 8 1 central region observed with the HERA receiver. Each spectrum has a velocity scale from -600 to 
600 km s _1 and a temperature scale from -20 to 70 mK. The positions are offsets relative to the M 8 1 nucleus assumed of coordinates 
RAy 2 ooo.o = 09' 1 55 m 33 s , DECy 2 ooo.o = 69° 03' 55". 



Table 3. R21 line ratio for spiral arms of M 81 by Brouille t et al.l 
(1991). 



Field 


Offsets 
(",") 


J r,„ A (i -o)dv 

[Kkms- 1 ] 


J T ml ,(2 - l)dv 
[Kkms- 1 ] 


A21 


Nl 


(0,0) 


1.14 


2.11 


1.85 


N5 


(0,0) 


1.52 


0.70 


0.46 




(10,0) 


1.27 


1.05 


0.83 




(-10,0) 


1.39 


0.70 


0.50 




(20,0) 


0.89 


0.53 


0.59 




(-20,0) 


1.52 


1.05 


0.69 




(20,-10) 


2.03 


1.05 


0.52 


N7 


(0,0) 


2.53 


1.23 


0.49 




(20,0) 


1.52 


1.93 


1.27 




(-20,0) 


2.28 


1.05 


0.46 




(0,20) 


0.51 


1.40 


2.77 




(0,-20) 


1.39 


4.04 


2.90 


S2 


(-10,0) 


0.76 


1.05 


1.39 




(0,-10) 


2.03 


2.28 


1.13 




(20,0) 


2.03 


1.40 


0.69 




(-20,0) 


1.65 


1.40 


0.85 



mal equilibrium at temperatures of 15 - 20 K. The J — 2 level 
of CO is populated at these temperatures, and Braine & Combes 
(1992) computed that the expected line ratios are between 0.7 
and 0.9. In addition, it is important to remember that the CO 
emission is not always thermalized and that the excitation tem- 
perature may change along the line of sight and within a single 
beam. 

We computed the R21 integrated intensity ratios for the off- 
sets observed in the M 81 center, after smoothing the 12 CO(2-l) 
data to the 12 CO(1-0) beam resolution (22") and correcting for 



the different beam efficiencies. In Table|2]the main results of the 
Gaussian fits of the observations are listed. In this table, columns 
(1) and (2) are the offsets (Aa and A/3) referred to the center of 
the galaxy, which coordinates are given in Table [T] columns (3) 
and (6) are the beam-average radiation temperatures (T m b) for 
the two transitions; columns (4) and (7) are the line-widths (Av) 
for the two transitions; column (9) is the ^21 = hxlha convolved 
line ratios; columns (5) and (8) are the CO intensities of the two 
transitions obtained from Gaussian fits. 

The first striking note is the low temperature both in CO(l- 
0) and CO(2-l) transitions, which characterizes the emission of 
all observed positions in the M 8 1 central region. For those de- 
tections > 2<t, the average brightness temperature in CO(1-0) is 
33.44 mK, while in CO(2-l) it is 21.47 mK. 

The /?2i ratios found in the central regions of M 81 are quite 
low: where the CO(1-0) emission is particularly weak, that of 
the CO(2-l) line appears weaker or difficult to detect, producing 
a low Ri\ value. The average 7?2i ratio for the center of M 81 
is 0.6 8, a value atypical of ga lactic nuclei. Actually, in CO sur- 
veys (Braine & Combes 1992) the central regions (inner kpc) of 
the galaxies have, on average, /?2i ~0.89, rar ely reaching higher 
ratios (e.g. for NGC 3310). According to iBraine & Combesl 
(1992), this happens because the inner molecular gas is opti- 
cally thick in 12 CO lines and is cool but not cold (T > 10 K). 
Our low ratio for M 81 is lower than the average value of 0.89 
for galaxies similar in morphology and distance. In addition, iso- 
lated galaxies appear to hav e line ratios that are low er than in- 
teracting/perturbed galaxies dBraine & Com bes 1992). M 81 ex- 
hibits a low value, despite being a classical example of interac- 
tion. 
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Integrated 12 CO(2-l) emission. The contour levels are be- 
0.410 to 4. 10 K km/s in steps of 0.410 K km/s. The beam 
is indicated at the bottom left. 
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Fig. 6. Integrated 12 CO(2-l) emission detected with the HERA 
receiver superposed on the 12 CO(2-l) emission found with A 
and B receivers. The contour levels are the same as shown in 
Fig. |5]between 0.410 to 4.10 K km/s in steps of 0.410 K km/s 
and the intensity levels detected with A and B receivers are in 
grey scale. The beam of 12" is indicated at the bottom left. The 
intensity peaks detected with the two observing configurations 
agree quite well between them. 



Studies of the disk in single galaxies show, in general, 
lower line ratios than those in the galactic nuc lei. IC 342 
dEckartet al.lll990l). NGC 6946 dCasoliet al.lll990l) . and M 51 
( Garcia-Burill o etalll 1 9931) are examples in which the ^21 ratio 
typically decreases from ~1 in the galactic nuclei to 0.4-0.6 in 
the spiral arms. 

To test if /?2i in M 81 shows a similar trend and the ratio 
decreases, which is necessary to use data in the literature. There 
is only one observational study on M 81 with data of the spira l 
arms in both CO(1-0) and CO(2-l) lines dBrouillet et al.lll991l) . 



We used Fig. 2 of iBrouillet et alj d!991l) to study only the fields 
with detections for the two lines and computed the correspond- 
ing /?2i ratios after the transformation from detected T* A to T m b, 
needed for the comparison with our data. The results of this anal- 
ysis are collected in Table[3] In this table, columns (1) and (2) are 
the fields a nd the correspo nding coordinates respectively defined 
in Fig. 2 of Brouill et et alJ d !991l) : columns (3) and (4) are the in- 
tensities l ines Zip and /21 respe ctively computed from the T* A de- 
tected by IBrouillet et all dl99ll) : column (5) is the R 2 i = /21//10 
line ratio. Table [3] shows that the R 2 \ line ratio spans in a wide 
range of values, without a clear radial trend. For instance, the 
regions N5 and N7, that are in the same area of the northeastern 
spiral arm, show very different values for R 2 \ ■ This confirms the 
"anomalous" R 2 \ behavior of M 81 with respect to other spiral 
galaxies. 

In conclusion, M 81 is a galaxy with low intensity and an 
anomalous behavior of CO lines. In this re spect, it is similar 
to the central CO-poor g alaxies NGC 7331 dYoung & Scoville l 
1982; lYoungetalj|l995l) . which has R 2X ~ 0.54 dlsrael & Baasl 
1999) and M 31, where similar low ratios ha ve been found for 
indivi dual dark clouds in the central region dAllen & Lequeuxl 
1993). In this case, the authors attributed this fact to a subther- 
mal excitation of the gas, at least for the CO(2-l) transition. 
These low ratios may be understood if coming from cool or cold, 
diffuse gas with a molecular hydrogen density (n(H2)) not high 
enough to excite the CO molecules. 



4. Clumping properties of the gas 

There is evidence that molecular clouds have a high degree of 
internal fragmentation on a ll sca les: this med i um is defined as 
clumpy. iRegan et al] d2001l) and iHelfer et all d2003l) . used the 
BIMA SONG survey -a systematic imaging survey in CO(1-0) 
emission of a sample of 44 nearby spiral galaxies- to study the 
molecular gas emission in the centers and inner disks of galaxies 
in their sample with a resolution of 6". Their CO maps display a 
remarkable variety of molecular gas morphologies, and the CO 
surface brightness distribution s hows substructures wi th dimen- 
sions of a few hundred parsecs. ISakamoto et al] d 19991) found a 
similar result with a CO aperture-synthesis of the central regions 
of 20 nearby spiral galaxies using the Nobeyama and Owens 
Valley millimetric arrays with a typical resolution of 4". In many 
galaxies they detected CO clumps of subkiloparsec size that they 
defined as giant molecular associations (GMAs). 

In Figs. |5] and |6] it appears difficult to distinguish individual 
structures, such as molecular clouds -probably present- and esti- 
mate their characteristic parameters (dimension, mass, etc.). To 
study the molecular gas structure on a small scale, it is necessary 
to use specific tools to identify single CO clouds. 

There are several different methods of decomposing clumps 
in a 3-dimensional space; i.e. position-position-velocity space, 
and all the techniques developed allow characterization of the 
clumpy structure: size, line width, and one of the most impor- 
tant parameters of a clump, the mass. The algorithms generally 
used to study the clumpy nature of the gas are GAUSSCLUMPS 
dStutzki & Guestedl990tlKramer et al.ll998l) and GAUSSFIND 
(Williams et all 11994 . GAUSSCLUMPS decomposes the 3- 
dimensional data cubes into a series of clumps, where each 
clump corresponds to a peak position and is assumed to have a 
Gaussian shape. GAUSSFIND works on the 3-dimensional data 
cubes as the eye would find a data set by using local maxima. 
Usually, the results of the two algorithms agree with each other. 
In general, GAUSSCLUMPS is able to more easily find low- 
mass clumps than GAUSSFIND; and since M 81 shows a weak 
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Number 


Aa 


A<5 


v 


To 
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Ax 2 
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Av 


dv/dr 


ip v 




["J 


["] 


[kms '] 


[K] 


["] 


["] 


[°] 


n -It 

[km s '] 


n -It 

[kms '] 


[°] 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


1 


6.9 


20.0 


6.13 


0.07 


15.8 


19.2 


83.7 


17.2 


1.67 


88.2 


2 


9.4 


8.8 


238.31 


0.06 


12.0 


12.1 


5.9 


12.1 


2.76 


-111.8 


3 


-14.8 


-16.1 


133.76 


0.05 


12.0 


12.1 


101.3 


12.1 


2.29 


-106.3 


4 


-7.1 


27.3 


132.12 


0.05 


12.0 


17.2 


94.5 


13.9 


2.05 


-69.7 


5 


4.4 


-20.2 


341.63 


0.05 


12.0 


12.1 


36.9 


12.1 


1.72 


-138.2 


6 


8.9 


2.7 


524.09 


0.05 


12.0 


19.0 


43.8 


14.3 


2.04 


-90.2 


7 


-17.3 


-20.6 


-21.47 


0.04 


12.0 


17.2 


48.5 


13.9 


1.40 


-86.3 


8 


9.2 


-28.2 


264.09 


0.04 


12.0 


12.1 


56.9 
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3.04 


-149.6 


n 

y 
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A A A 
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10 


-14.7 


-27.0 


263.61 


0.04 


12.0 


12.1 
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12.1 


4.41 


28.6 


n 


7.9 


9.0 


30.24 


0.04 


12.0 


13.0 


121.2 


12.5 


3.42 


-33.7 


12 


1.6 


27.2 


56.42 


0.04 


12.0 


23.6 


88.8 


15.1 


1.39 


-105.3 


13 


8.4 


-15.5 


316.24 


0.04 


12.0 


12.1 


87.7 


12.1 


5.83 


4.7 


14 


9.2 


-14.9 


82.50 


0.04 


12.0 


12.1 


22.6 


12.1 


3.62 


115.8 


15 


-32.9 


-9.2 


-22.05 


0.04 


12.9 


20.0 


81.0 


15.3 


1.04 


-18.8 


16 


-14.5 


21.9 


-22.41 


0.04 


19.0 


12.1 


29.3 


14.5 


2.18 


42.0 



CO emission, we decided to use GAUSSCLUMPS. In addition, 
this algorithm can work with the package GILDAS, already used 
here. 



4. 1 . Identification of molecular associations 

The first application of GAUSSCLUMPS to CO(2-l) data ob- 
tained using HERA receiver found 27 molecular associations in 
the central 1.3 kpc; however, it is necessary to fix the constraints 
to check misidentification. The algorithm forces the unresolved 
clumps to have a size corresponding to the spatial and veloc- 
ity resolutions, so if the 3D size of a clump is exactly equal 
to the resolution in each axis, this clump could be an artefact. 
With a beam of 12", only clumps with the intrinsic FHWM size 
A£ > 10" are realistic, so all the other clumps have been dis- 
carded. The second constraint is on the intrinsic brightness tem- 
perature, which must be at least 5 times the noise l evel of the 
original map, T > 0.035 K (IStutzki & Guestenlll990l) . After ap- 
plying these constraints, the number of molecular associations 
found in the M 81 central region has become 16 and their prop- 
erties are listed in Table [4] In this table, column (1) is the GMA 
number; columns (2) and (3) are the GMA center positions re- 
ferred to the center of the galaxy, which coordinates are given in 
Table [TJ column (4) is the GMA center velocity; column (5) is 
the GMA peak intensity; column (6) and (7) are the GMA spatial 
FWHMs along the two principal axis; column (8) is the orienta- 
tion of the first principal axes; column (9) is the GMA velocity 
FWHM; column (10) is the GMA internal velocity gradient; col- 
umn (1 1) is the direction of internal velocity gradient. 

Considering that the Gaussian-shaped molecular associa- 
tions have two principal axes, we find a mean intrinsic di- 
ameter for all molecular structures identified of 14" (~250 
pc). Molecular structures with these dimensions are defined as 
GMAs. The maximum GMA peak intensity is 0.070 K, the mean 
intensity is 0.046 K and the mean velocity FWHM 13.5 km s _I . 

Summing up the area corresponding to the intrinsic FWHM 
extent of the fitted GMAs in the selected field, one obtains a total 
area covered by GMAs that is from 2.2 to 4.5 times smaller than 
the field, therefore the filling factor of the GMAs is low, about 
0.3. 



4.2. GMAs masses 

For an optically thin gas, or at the other extre me for a com- 
pletely optically thick gas dKramer et alJll998h . the mass of a 
molecular association is proportional to the fitting parameters 
of GAUSSCLUMPS: the molecular association intensity, the 
molecular association velocity width, and the spatial FWHMs 
along both two principal axes: 



M oc T pk x Av x Axi x Ax2. 



(1) 



We now present mass computations in the two hypotheses 
below. For the estimation in the optically thin case, we adopt the 
Rayleigh- Jeans approximation, where: f(T ex ) = Y=^^tjW^\ ~ 
T ex , since T ex > 15 K. LTE is also assumed, meaning that 
T ex is the same for all levels. The total CO column density 
Nco is related to the population Ni of level i by Ni/Nco = 
giexp [-Ei/lcT ex ] IQ, where g, and Ej are the degeneracy fac- 
tor and the energy of the level i, and Q the partition function is 
Q(T) = 2kT/hvi- = T(K)I2.16. This total CO column den- 
sity, derived from the integrated intensity of the CO(2-l) line, is 
given by the following equation, in the optically thin limit: 



N, 



{mol/ cm 2 ) 



1.2xl0 l3 xr„W I T mb (K)dv(bn/s). (2) 



The corresponding molecular mass (taking into account helium) 



M 1.2x10^,^x7^) 



(M ) 



(M ) N( l2 CO) 
J j T mb (K) dv(km/s)dA(cm 2 ). 



(3) 



In terms of fitting parameters in Table [4] (see IChi & P ark 2006) 
and when assuming an excitation temperature of T ex 15 K, an 
abundance ratio N(H 2 )/N( 1Z C0) - 10 4 , and a distance of 3.6 
Mpc, Eq. (0, valid in the optically thin case, becomes 



I hin 



T mb Av Axi Ax 2 , _ . 
x x — — x — — x 10.5. 



(M Q ) (K) (km/s) (") (") 



(4) 



In this hypothesis, the ratio between the observed peak bright- 
ness temperature and the excitation temperature is the filling 
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Table 5. Properties of the GMAs derived by application of GAUSSCLUMPS. 
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10 3 


10 5 


10 6 




10 21 




[pc] 


[M B ] 


[M B ] 


[M ] 




[mol cm" 2 (K km s~ 1 )~ 1 ] 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


1 


151.88 


3.83 


7.20 


8.54 


11.86 


2.29 


2 


105.08 


1.11 


2.08 


2.92 


14.07 


2.71 


3 


105.08 


0.92 


1.73 


2.92 


16.89 


3.26 


4 


125.28 


1.50 


2.83 


4.60 


16.27 


3.14 


5 


105.08 


0.92 


1.73 


2.92 


16.89 


3.26 


6 


131.67 


1.71 


3.21 


5.12 


15.93 


3.07 


7 


125.28 


1.20 


2.26 


4.60 


20.34 


3.92 


8 


105.08 


0.74 


1.38 


2.92 


21.11 


4.07 


9 


123.81 


1.17 


2.19 


4.48 


20.43 


3.94 


10 


105.08 


0.74 


1.38 


2.92 


21.11 


4.07 


11 


108.91 


0.82 


1.54 


3.23 


21.04 


4.06 


12 


146.74 


1.79 


3.37 


6.36 


18.86 


3.64 


13 


105.08 


0.74 


1.38 


2.92 


21.11 


4.07 


14 


105.08 


0.74 


1.38 


2.92 


21.11 


4.07 


15 


140.06 


1.66 


3.11 


6.23 


20.03 


3.86 


16 


132.22 


1.40 


2.63 


5.28 


20.10 


3.88 



factor of the porous molecular associations, on the order of 3 
xlO 3 . Applying Eq. |@}, we find very small molecular associa- 
tion masses contained in a low mass ranee, from 7.37 x 1O 2 M 
to3.83xl0 3 M o , with a mean value of 1 .3 1 x 1O 3 M . The masses 
computed for all molecular associations are listed in Table [5] In 
this table, column (1) is the GMA number as defined in Table 
IH column (2) is the effective radius of the GMAs; columns (3), 
(4), and (5) are the masses of the GMAs computed in the opti- 
cally thin limit, in the optically thick limit, and using the virial 
theorem, respectively; column (6) is the ratio a,.,,. = M V j r /M b s 
for each GMA, where M t, s = M t hi C k', column (7) is the derived 
X conversion factor for each GMA. 

In the optically thick case, the molecular mass of the molec- 
ular associations is determined by the H2 density (Nh 2 ) using the 
standard X conversion factor value (X = 2.3 x 10 20 mol crrT 2 (K 
km s" 1 )" 1 . IStrong et al.ll 19881) . N Hl = X JY CO (i-o)dv. Since we 
resolve single GMAs in the CO(2-l) transition, using the mean 
^21 = hxlho value computed here (7?2i = 0.68), the H2 density 
can be expressed by 



x 1 r 

ollcnf_\ R 2l J 



Tcoa-i) (K) dv (km Is) 



(mol I Cm 2 ) ( mot/cm 2 \ 

v \ Kkm/s I 

= 3.4 x 10 20 J" T C0 (2-i) (K) dv (km/s) , (5) 
and the corresponding molecular mass is 
M -3.4xl0 2 ° m " 2 



(M ) 



(M ) 



x J J Tco(2-x)(K)dv(kmls)dA(cm 2 ). (6) 

In terms of fitting parameters in Tableland assuming a distance 
of 3.6 Mpc, Eq. (0 gives the GMA mass in the optically thick 
limit: 



M thick T mb Av Axi Ax 



• x 



(M ) (K) (km/s) (") (") 



x — - x —4 x 1.97 x 10 3 



(7) 



Applying Eq. (0, the GMA masses range from 1.38 x 10 5 M G 
to 7.20 x 1O 5 M , with a mean value of 2.46 x 10 5 M o . The total 



mass of the GMAs is 3.22 x 1O 5 M . The masses computed for 
all GMAs are listed in Table g] 

That the masses obtained through the optical thin hypothesis 
are much lower than those calculated in the optically thick case 
is not "anomalous": instead, it is expected that the optical thin 
hypothesis underestimates the masses, since GMAs have a rich 
substructure of subclumps with high column density. The sub- 
clumps are probably optically thick, as is generally the case in 
external galaxies, as suggested by the high 13 CO emission com- 
monly observed. Anyway we decided to also analyze the opti- 
cally thin case because M 8 1 represents an extraordinary labora- 
tory where all possibilities should be considered even if unusual 
for similar galaxies. 

4.3. Virial equilibrium 

An alternative way of determining GMA masses is to compute 
their virial mass. Assuming a balance between the total kinetic 
energy T of a GMA and the gravitational potential £2, and ne- 
glecting external pressure, the virial mass M V! > of the clouds can 
be calculated using the formula 



KxDrX Av 2 



(8) 



where M V! > is in M Q , D e is the effective diameter in pc of the 
equi valent sphere of radiu s R = (AxiAjc2/4) 05 , and Av is in km 
s- 1 (Solomon et al. 1987). We assume that the cloud is a sphere 
with a power-law (a — I) density distribution. In this case, the 
constant K in formula © is equal to 95. The rougher assump- 
tion of uniform density © would have instead produced K= 104. 
Here, the hypothesis of spherical shape for the cloud makes the 
formula independent of the galaxy inclination. 

A comparison between the virial mass and the mass com- 
puted from the observed intensity lines (M b s ) in the optical thick 
hypothesis is informative for the physical state of the clouds, 
their optical thickness, and for the actual X= N(H2)/I(CO) con- 
version factor. 

The virial masses of the 16 GMAs we identified span from 
2.92 x 10 6 to 8.54 x 1O 6 M , with a mean value of 4.31 x 1O 6 M . 
In calculating M V! >, we have neglected the mass of the atomic 
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hydrogen (M H1 = 0) since lAllen et al.1 d!997l) . studying the HI 
profile in M 81, found a deficiency of atomic hydrogen in the 
central ~3 kpc, and the HI emission only begins to be significant 
in the inner spiral arms beyond ~5 kpc from the nucleus. 

We also checked that the stars do not contribute significantly 
to the mass inside the volume of a GMA. In the hypothesis that 
the H2 mass is computed from the virial mass, resolved GMAs 
have mean surface densities ranging from 84 M /pc 2 to 118 
M /pc 2 with a mean value of 91 M /pc 2 . The mean value of 
the mass-to-light ratio of visible m atter expected for M 81 is 
M/L B = 5.5 M /L (see Te niesll 19941) . From the exponential sur- 
face brightness distribution of the stellar disk of the galaxy, we 
then deduced a mean stellar surface density of 36 M /pc 2 in the 
region. In addition, if we consider that GMAs normally have a 
thickness of about 50 pc while the stars are distributed in regions 
that are at least ten times thicker, the actual mean volumic stellar 
density is negligible with respect to the molecular density inside 
a GMA. 

All the GMAs have a virial mass » of the mass computed 
from CO observations, and the mean ratio a,,,> = M V j r /M i, s com- 
puted for our GMAs is a,,,> ^ (10 - 20) » 1. The a v , r ratios 
computed for all the GMAs are listed in Table|5] 

4.4. TheN(H 2 )/Ico conversion factor 

The value of the conversion factor between the Hi column den- 
sity of the molecular hydrogen NiHi) and the CO integrated 
intensity (X = N{H2)I Ico) is still subject to debate. Usually 
the standard value of X found for the Milky Way is universally 
adopted for all the galaxies and for all the regions of the same 
galax y (X = 2.3 x 10 20 mol cm" 2 (K km s 1 ) 1 . IStrong et al.l 
1988). This assumption often turned out not to be completely 
correct. 

The X conversion factor is determined by various factors, 
such as the metallicity, the temperat ure, the cosmic ray den- 
sity, and the UV r adiation field (see iMalonev & Blackl fl988t 
Boseljietalj|2002). For instance, galaxies with lower metal- 
licity generally have a higher X factor; when studying 5 local 
group galaxies, IWilsonl d 1995b found that the conversion factor 
increases by a factor of 4.6 when metallicity decreases of a factor 
10. The X value also changes with the morphological type of the 
galaxies: galaxies earlier than Scd show values comparable to, 
or lower than, the Galactic one, while extremely l ate-type spirals 
or irr egular galaxies tend to show higher values (Nak ai & Kunol 
1 1995b . A dependency of the X conversion factor on Ico is also 
suspected: galaxies or observed positions with low Ico (Ico < 
20-40 K km s _1 ) show a higher X. An example of this behavior 
is the galaxy M 51, where Ico < 20 K km s" 1 and the X derived 
valu e is a factor 2.6 lower than the mean Galactic value (see Fig. 
3 of iNakai&Kunoll 19951) . 

Since M 8 1 shows atypical properties for its molecular gas, 
the CO to H2 conversion ratio is also suspected of departing sig- 
nificantly from the mean Galactic value and in general from that 
of galaxies with the same morphological type. Thus, M 8 1 rep- 
resents an exceptional laboratory to explore the variation of the 
X-ratio. 

We derived the X conversion factor for all individually re- 
solved molecular associations assuming the virial equilibrium 
for single molecular complexes and using the large velocity gra- 
dient (LVG) approximation for the format ion of the CO line 
dYoung & Scovilldll984l:ISanders et al.lll984l) . This method has 
also been used to derive the conversion factor in nearby spi- 
ral galaxies such as M 31, M 33, and the Magella nic Clouds 
dWilson & Scovilld[l 990. 19921: lArimoto"iraT][l996h . 



In practice, we computed the X ratio (=N(H2)/Ico) deriving 
the H2 density from the virial masses (N(H2) = M V j r /(niH 2 xdA), 
where is the mass of the hydrogen molecule and dA the area 
occupied by the cloud), and Ico is the CO observed intensity. For 
the individually resolved GMAs of M 81, assumed virialized, 
we found that the X conversion factor takes a much higher value 
than the mean Galactic value for all clouds, ranging from ~10 
to ~18 times higher than the standard X and a mean value of 
X = 3.6 x 10 21 mol cirT 2 (K km s -1 ) -1 . The X conversion factor 
values computed for all the GMAs are listed in Table [5] 

Our X factor values are high, not only if compared to the 
standard Galactic factor, but also with respect to values found 
in other nearby galaxies where the X ratio has been obtained 
with the same method we used. Studying a sample o f nearby 
spiral and dwarf irregular galaxies, lArimoto et alj d 19961) derived 
X values using CO data found in the literature and computed 
the corresponding virial mass. They found galaxies with both X 
values lower and higher than in the Galaxy, such as M 3 1 with 
X = 1.8 x 10 20 mol cm" 2 (Kkm s" 1 )" 1 , M 33 with X = 4.1 X 10 20 
mol cm" 2 (K kms" 1 )" 1 , and IC 10 where X = 6.6 x 10 20 mol 
crrT 2 (K km s Even if they found high X values, these H2- 
to-CO co nversion factors are s maller than our X ratios obtained 
for M 8 1 . lArimoto etail d 19961) found a X ratio (=3. 1 x 10 21 mol 
cirT 2 (K km s -1 )-') only for the SMC within the same order 
of magnitude of M 81. They specify that this value has been 
obtained by considering virial masses and CO luminosities for 
larger-scale complexes in the SMC and that this high X would 
be due to a smaller amount of CO in a diffuse medium that is 
metal dependent, as well as photo-dissociated by the strong UV 
field, and can be present only in dense molecular clouds. 



5. Discussion 

5. 1 . Comparison with previous works 

Our molecular gas detections in the center of M 8 1 agree with the 
only two works present in the literatur e that report CO emission 
in the nuclear region of this galaxy (Sage & West pfahl 11991b 
Saka moto et afl l2001l) . In comparison with th e giant molecular 
clouds (GMCs) in the oute r regions of M 81 (Tayl or & Wilson! 
119981: IBrouillet et~ai1ll998l) . the CO emission that we detect in 
the M 81 center is particularly cold and unusual. The masses 
we find for GMAs in the CO(2-l) map are similar t o the higher 
masses observed in M 33 dWilson & Scovill e 1990), and higher 
than what is measured in the SMC flRubio et al.lll993l) . 

Our first result is that the nuclear region of M 8 1 (within 
~300 pc) appears de void of CO emission, in a greement with 
IBrouillet et a Ll d 19881) and lSakamoto et ail fcOQl). However, we 
detect CO em ission outside the 300 p c radius. This supports 
the claims of ISage & Westpfahl d!991l) who, found CO emis- 
sion when observing one offset towards the center galaxy with 
the NRAO 12m telesco pe and a beam of 55" (see Fig. 1 of 
ISage& Westofahlll991l) . There is a good correspondence be- 
tween our velocity profile and theirs. 

A more useful and interes ting comparison can b e done be- 
tween our results and those of ISakamoto et al.l d200ll) . In detail, 
the "pseudoring" that we detected at ~ 460 pc from the nucleus 
in the north east direction (F ig. [3]) corresponds to the maximum 
intensity of Sakamoto et al.1 (bOQll) . Obser ving the central kilo- 
parsec with the Nobeyama Millimetric Array-NMA (resolution 
of 6.9 "x 5.8 " that corresponds to ~ 120x100 pc at the distance 
of M 81), they produced maps only for the 12 CO(1-0) line. Since 
they made interferometric measurements, with an FWHM pri- 
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Table 6. Comparison of the results obtained applying the aperture photometry at the FUV-GALEX image of M 8 1 in different 
regions observed in CO(l-O). 



Reference 


Ico 


Ffuv 


Go 


IcofFpuv 






io- 14 




10 14 




[Kkms- 1 ] 


[erg cirr 2 s _1 A -1 ] 




[K km s"V(erg cm- 2 s" 1 A" 1 )] 


Brouillet etal. (1998) 


25.90 


1.37 


0.51 


18.88 


Knapen et al. (2006) 


20.26 


1.48 


0.35 


13.65 


This paper 


53.50 


1.13 


0.31 


47.18 



mary beam of 1' and a beam of ~ 7", extended emission seen 
with the 22" IRAM CO(1-0) beam is filtered out. 

In spite of this, the total mass of molecular gas that we 
find, a ssuming the same standard X-fac tor value of [Strong et al] 
(1988), is quite comparable to that of Sakamo to et al.l (1200 ll) : 



our resolved GMAs in CO(2-l) emission have a total mass of 
3.22 x 10 6 M B , only a l ittle bit lower than ~ 1 xlO 7 M found by 
Sakamoto et al.l (12001 1) for the "p seudoring" in CO(l-O) emis- 
sion. In addition to the result of Sakamoto et al. (2001), who 
found the molecular gas mainly located in a lopsided distribu- 
tion in the northeast direction, we also detected some emission 
towards the southwest. Our peak temperatures for the I2 CO(l- 
0) line are very low, ranging from ~15 to ~64 mK, lower than 
the avera ge peak brightness te mperature (~0.3 K) of each GMA 
found by Sakamoto et al. (2001). This is compatible with the dif- 
ferent beam sizes, the synthesized beam of the interferometer 
being 10 times smaller in surface. 

Another useful comparison can be made with molecular 
clouds identifications in other external galaxies. The masses 
(~ 10 5 M o ) and the diameters (~ 250 pc) of the CO clouds iden- 
tified in M 81 are not much unu sual for giant galaxies of the 
same type dTavlor & Wilsonll998l) . Even if in general the GMCs 
and GMAs in the galaxies of the Local Group are very similar, 
there is evidence of variations in the properties of the molecular 
clouds of galaxies of different morphologies, which range from 
early-type spirals (M 31 and M 81) to late-type spirals (M 33 
of Scd morphological type) and from dwarf irregulars (IC 10) 
to dwarf ellipticals (NGC 205). For the variety of characteristics 
found in clouds of gala xies of these sub-morphological types, 
iTavlor & Wilsonl £1998) conclude that in general very late-type 
galaxies should have smaller and less massive GMCs than early- 
type galaxies. Our results agree with this conclusion. 

Finally, we know that, in the Local Group and in sys- 
tems in interaction with the Local Group, macroscopic prop- 
erties of molec ular clouds are related by power relationships, 
first noted by iLarsonl (1 1 98 lh an d then studied in det ail by 
ISolomon etal] dl987l) .In particular. lSolomon et all d!987l) found 
a tight relationship between the Galactic dynamical molecular 
cloud masses, measured by the virial theorem, and the CO lu- 
minosity (M VJ > oc (Leo) ' 81 )- For several external galaxies, this 
power law is even valid if with little variations of the power 
index dRand & Kulkarnilll990HRandlll995l; lRosolowskv & Blitz 
120051) . It is interesting to test if M 81 also shows non-typical 
behavior in this case. Our very weak CO detections correspond 
to CO luminosities of the order of magnitude 10 4 K km s" 1 pc 2 
with virial masses within an order of magnitude of 10 6 M Q (see 
Table |5). Comparing our results with the Fig. 2 of lSolomon et al] 

~~ 53>T' 



dl987l) 7we can see that with the Leo — 10 4 K km s" 1 pc 2 virial 
masses we infer are overestimated by more than one order of 
magnitude with respect to the typical trend in our Galaxy, but 
also in external galaxies. This result agrees with the indications 
given by the ar vi > ratios calculated in Sect. 4.3, and in general it 



confirms the unusual scenario in the molecular gas properties of 
M81. 



5.2. The X-factor problem 

The X-factor peculiarity does not appear to be due to a particu- 
lar metallicity, since M 81 has a metallicity of Z=0.03, without 
any evident gradi ent from the central region to the outer disk 
dKong et al.ll2000l) . This "normal" metallicity value would pro- 
duce a X conversion factor similar to that of the Milky Way. 
Some regions in M 81, however, have higher metallicity and are 
mostly located in the spiral arms and HII regions. But in general, 
our estimate of the X-factor appears significantly high. 

A possible source of error may be the assumption of a fully 
virial equilibrium for the observed clouds. Under this hypothe- 
sis, we calculated the X-factor using X = M vlr jLco, where the 
M V i r has been computed with Eq. ©, and Leo is the CO lumi- 
nosity derived from our observations. Application of the virial 
theorem to complex molecular structures may not be completely 
appropriate, even if the expected deviation cannot be higher than 
a factor ~2. However, if the GMAs are not virialized and the 
standard X factor is applied, the molecular content of M 81 
would be exceptionally low, so M 8 1 would be peculiar in both 
cases. 



5.3. Heating of the gas 

The reason for the non-typical results on the molecular gas in 
M 8 1 is probably the exceptionally low CO excitation temper- 
ature that seems to characterize the center of this galaxy. Why 
is the molecular gas absent or subthermally excited in the M 8 1 
center? 

A possible link could exist between low CO emission and 
the lack of ultraviolet emission that appears to characterize the 
M 8 1 central region. The ultraviolet image of M 8 1 realized with 
the satellite Galaxy Evolution Explorer (GALEX) clearly shows 
a lack of recent star formation, traced by the UV emission, in the 
nuclear region that we observed in CO. 

The existence of a link between molecular gas and far- 
ultraviolet (FUV) emission can be justified considering that the 
interstellar medium is excited, dissociated, and ionized by the 
FUV photons produced by young O and B stars. Atomic gas 
in the ISM recombines into molecular form mainly through the 
catalytic action of dust grain surfaces. Hydrogen nuclei cycle re- 
peatedly from the molecular (H2) to the atomic phase and back 
again at rates that depend on the incident FUV flux, the total 
volume density of the gas, and the dust-to-gas ratio. The action 
of the FUV flux that interacts with the surface of the molecular 
clouds can strongly affect their physics and chemistry by heat- 
ing, dissociating, and exciting the gas. 

In regions of the ISM where the physics is dominated by 
FUV photons -called photodissociation regions (PDRs)- useful 
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information can be obtained on the physical state of the gas from 
observations of the spectral lines emitted by the excited atoms 
and molecules in those regions. 

In M 81 the possibility that low CO emission a nd weak FUV 
flux ar e correlated has already been proposed by Knapen et al. 
(120061) . When observing a region of the western spiral arm in 
M 8 1 , they explained the lack of or very weak CO( 1 -0) emission, 
observed far from FUV sources, as a consequence of insufficient 
heating and excitation of the molecular gas. Usually, a lack in the 
CO emission is instead interpreted as due to the absence of mea- 
surable quantities of molecular gas or due to a low density gas 
not heated sufficiently to detect it. The limit of this interpreta- 
tion, based on the linear proportionality between the CO surface 
brightness and the column density of the molecular hydrogen, is 
to neglect ef fects linked to local h eating mechanisms. If the ex- 
planation of Knap enet al.l d2006l) is correct, the flux of cosmic 
rays and the FUV surface brightness would be too low to heat 
the ISM, and in particular, the molecular component. It is pos- 
sible that this correlation of low-FUV/low-CO is also applicable 
to the central region that we observed. 

To verify this interpretation, we performed aperture photom- 
etry in the FUV emission on three fields already observed in 
CO: our c entral region o f 1 .6 k pc, the field of spiral arms ob- 
served by Brouille t et al. ( 19981). and th at - again on the spiral 
arms - observed by Knapen et al. (200o]). The fields selected o n 
the spiral arms are the same of Fig. 2 of lKnapen etaTJ ([2006). 
We used the FUV-GALEX image of M 8 1 expressed in intensity 
and already sky-subtracted and performed the aperture photom- 
etry using the task QPHOT of IRAF, able to do a quick aper- 
ture photometry. This procedure allowed us to know the FUV 
flux emitted in regions observed in CO(1-0) and hence to com- 
pute the ratio between the two emissions in different fields, but 
of similar dimensions, of the galaxy. The aperture photometry 
was done with a rad ius of 40" for our c entral region and for 
the one observed by Knape net al.l d2006l) with the Nobeyama 
Radio Obser vatory (NRO), and wi th a radius of 32" for the field 
observed by Broui llet et alj (1998). The results of the aperture 
photometry are summarized in Table|6] In this table, column (1) 
gives the references relatively to CO(1-0) emission; column (2) 
is the CO(1-0) integrated intensity; column (3) is the FUV flux 
obtained by the aperture photometry; column (4) is FUV flux 
expressed in terms of the Go Hading unit; column (5) is the ra- 
tio between the CO(1-0) integrated intensity and the FUV flux. 
We express the FUV fluxes in terms of the Hadin g unit, since it is 
one of the units usually used to model PDRs (e.g. Kaufma n et alj 
1999; lAllen et alj|20 04). In these models each one-dimensional 
semi-infinite slab of gas with constant density of H nuclei is sub- 
jected to an equivalent one-dimensional flux of FUV photons Go 
measured in units of 1.6 x 10~ 3 ergs cirT 2 s over the photon 
energy range 6 - 13.6 eV. In these unit s, the average interstellar 
radiation field over An sr is Go = 1.7 (Draind ll978HAllen etalJ 
120041) . 

The central region emits a weak FUV flux, within the same 
order of magnitude as that emitted in the two fields of the spiral 
arms, and the higher ratio between CO and FUV emissions in 
nuclear region is mainly due more to the higher molecular gas 
intensity. If the low FUV emission is not able to heat the molec- 
ular gas sufficiently in the spiral arms, in the central region at 
the same quantity of FUV flux corresponds a molecular gas dis- 
tributed in substructures, but subexcited. Probably the central re- 
gion shows a higher CO/FUV ratio because in the spiral arms the 
molecular gas is less excited than the nucleus. 



6. Conclusions 

The physics of the molecular gas suggested by the observations 
reported in this paper shows M 8 1 to be a galaxy that is not only 
CO-poor in its central region, but that also has molecular clouds 
with unusual properties, if compared to galaxies with similar dis- 
tances and morphological type. The absent or very weak molec- 
ular gas emission in the nuclear region, the low 7?2i line ratio, and 
a particularly high X conversion factor value are its main pecu- 
liarities. The N(H2)/Ico conversion factor has a mean value of 
X = 3.6 X IP 2 1 mol cm" 2 (K km s" 1 )" 1 , 15.6 times higher than 
the X value of lStrong et al.l JT988) derived for our Galaxy. 

All unusual results found for the molecular component in 
the M 81 center are probably due to an excitation process of 
the gas non-typical for giant spirals in the Local Universe. The 
spiral arms of this galaxy are CO-poor, and it has been suggested 
that the FUV surface brightness is too low to heat the molecular 
gas component. The low CO emission we found in the center 
suggests that the gas in the nucleus is also subexcited, with even 
weaker FUV emission. We conclude that the lack of excitation 
of the gas, more than the absence of molecular gas, is the cause 
of the low CO emission there. 
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